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Using first principles density functional calculations, together with exact diagonalization of Fe- 
Mo Hamiltonian constructed in a first principles Wannier function basis, we studied the electronic 
structure of La doped double perovskite compound Sr2FeMo06- Our calculation show stabilization 
of kinetic energy driven antiferromagnetic phase for La rich compounds, in agreement with the 
results obtained on the basis of previous model calculations. 
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I. INTRODUCTION 

Double perovskites with a general formula A 2 BB'C>6 
where B and B' are transition metal ions and A is a 
rare-earth or alkaline-earth ion, are materials that have 
attracted enormous amount of attention in recent time 
due to the diversity of their applications as for ex- 
ample in the field of spintronics (Sr 2 FeMoOjS), mul- 
tiferroicity (Ba 2 NiMnOj2l), magnetodielectric materials 
(La 2 NiMnO E pH), magneto-optic devices (Sr 2 CrOs0 6 , 
Sr 2 CrReOj^). The choice of B and B' ions, provide the 
tunability of B-O-B' interaction, giving rise to a vari- 
ety of magnetic properties like ferromagnetism, antifer- 
romagnetism, ferrimagnetism and electronic properties 
like metallic, half-metallic and insulating^J. The pres- 
ence of two transition metal ions instead of one as in 
perovskite material is expected to give rise to far more 
tunability and richness of properties compared to simple 
perovskites. 

Perhaps the most studied member of this series that 
arose much interest is Sr 2 FeMo06 (SFMO). This mate- 
rial was report exhibit a large magnetoresis- 
tance(MR) effect with a fairly high ferromagnetic transi- 
tion temperature of about 410 K, opening up the pos- 
sibility of designing spintronics materials operating at 
room temperature. However, unlike colossal magne- 
toresistive compounds as manganites, this MR does not 
arise from electron-phonon interactions. Rather, it is 
extrinsic, of tunnelling magnetoresistive (TMR) origin. 
Since the report of the large MR effect and high mag- 
netic transition temperature, a number of experimental 
studies like NMRpl^CES^S, Hall measurement^, mag- 
netic measurement*^ have been carried out to charac- 
terize various properties of this material. There have 
been also a number of theoretical studies involving 
both first-principles calculation^U-Hfil as we h as model 
calculationg^IHSS] The unusually high ferromagnetic 
transition temperature in Sr 2 FeMo06 and related ma- 
terial like Sr 2 FeRe06 was rationalized 17 ^ 23 ! m terms of a 
kinetic energy driven mechanism which produces a nega- 
tive spin polarization at otherwise nonmagnetic site like 
Mo or Re. Following this idea, a double-exchange like two 
sublattice model was introduced and studied by differ- 



ent groupjgE2 1 S P B p gi while most of the studie^UmHH 
were restricted only to ferromagnetic phase, some of the 
studies^^H were extended to other competing magnetic 
phases too. Very recently^, the problem has been stud- 
ied in detail in terms of a full numerical solution of spin- 
fermion model and as well as in terms of reduced, clas- 
sical magnetic model. These studies predict that when 
the competing magnetic phases are taken into account, 
the electron doped model systems beyond a certain dop- 
ing prefers to have antiferromagnetic(AFM) arrangement 
of Fe spins compared to ferromagnetic(FM) arrangement 
of the undoped system. The predicted antiferromagnetic 
phase in electron-doped system is kinetic- energy driven 
rather than super-exchange driven, as is the case for 
example in Sr^FeWO^^, which is an insulating anti- 
ferromagnet with Neel temperature of sa 20 K. The su- 
perexchange driven antiferromagnetic phase is necessar- 
ily insulating while the kinetic energy driven AFM phase 
may not be so. The prediction of such an antiferromag- 
netic phase of different origin is therefore of significance. 
While the kinetic energy driven antiferromagnetic phases 
have been suggested in hole-doped rare-earth manganites 
(eg.the CE phase at half-doping^ , to the best of our 
knowledge, till date no reports of such analogous phases 
in double perovskites exist, thereby, opening up the pos- 
sibility of experimental exploration in this front. How- 
ever the afore-mentioned model calculations were carried 
out in two dimension and with single band, which was 
justified by the assumption that the dominant nearest- 
neighbor B-B' interactions are operative between orbitals 
of same symmetry and within a given plane. These re- 
strictions are not strictly true. Furthermore, the mag- 
netic ordering in real material is three-dimensional. A full 
three-dimensional, all orbital calculation without these 
approximations, is therefore necessary to put the possi- 
ble existence of the AFM phase in firm footing. 

Considering the above mentioned points, it is therefore, 
of interest to study the problem of electron doping using 
first-principles, density functional theory (DFT) based 
calculations. The first principles calculations which take 
into account all the structural and chemical aspects cor- 
rectly is expected to provide more realistic scenario and 
verification of predictions made by model calculations. 



The Sr ions in SFMO can be substituted for trivalent 
cations, like La, leading to S^-^LazFeMoOg- This would 
cause electron doping in the system, with 1 + x electron 
per formula unit in the conduction band, compared to 
1 electron per formula unit in the undoped SFMO sit- 
uation. To our knowledge, there exists very few first- 
principles study of the La-doped SFMO system. Few 
studies^SUlil that exist explored only the ferromagnetic 
phase, did not consider the other competing magnetic 
phases and were restricted mostly to Sr-rich part of the 
phase diagram. Motivated by the findings of the model 
calculations^, we considered it worthwhile to span the 
whole concentration range from x — 0.0 ie S^FeMoOg to 
x = 2.0 ie La2FeMo06 and study the relative stability of 
the various magnetic phases as one increases the carrier 
concentration through the increased doping of La. 

We have carried out our study both in terms of full 
ab-initio calculations as well as in terms of solutions of 
multi-orbital, low-energy Hamiltonians defined in a first- 
principles derived Wannier function basis. The structural 
optimization and total energy calculations of various 
magnetic phases have been carried out using the plane 
wave pseudopotential method as implemented in the Vi- 
enna Ab-initio Simulation Package (VASP)^, while the 
doping effect in first-principles calculations has been sim- 
ulated through supercell technique. The construction of 
low-energy Hamiltonian in first-principles derived Wan- 
nier function basis has been achieved through muffin- 
tin orbital (MTO) based Nth order MTO (NMTO)- 
downfolding technique^. The constructed multi-orbital, 
spin-fermion Hamiltonian defined in the first-principles 
derived Wannier function basis has been solved by means 
of real space based exact diagonalization technique. 

The rest of the paper is organized in the following 
manner. Section II contains the details of the employed 
methods and calculations. Section III is devoted to re- 
sults which consist of three subsections: (A) Total energy 
calculations, electronic structure and relative stability of 
various magnetic phases in doped compounds (B) De- 
termination of low-energy, few orbital Hamiltonian by 
NMTO-downfolding (C) Calculations of magnetic phase 
diagram and magnetic transition temperatures in terms 
of low-energy Hamiltonian. The paper concludes with 
section IV containing discussion and summary. 



II. METHODS AND COMPUTATIONAL 
DETAILS 



The first-principles DFT calculations were carried 
out using the plane wave pseudopotential method 
implemented within VASP. We considered exchange- 
correlation functionals within generalized gradient ap- 
proximation (GGA)P1 and GGA+lj2^. We used projec- 
tor augmented wave (PAW) potentials^ and the wave- 
functions were expanded in the plane wave basis with a 
kinetic energy cut-off of 450 eV. Reciprocal space inte- 
gration was carried out with a k-space mesh of 6 x 6 x6. 
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Figure 1: The Fe sublattice ordering of S^-aTa^FeMoOe. 
Shown are the A-type(left panel) and G-type (right panel) an- 
tiferromagnetic arrangement of Fe spins. In case of A- type an- 
tiferromagnetism the Fe spins in-plane are ferromagnetically 
coupled, while Fe spins between two adjacent plans are an- 
tiferromagnetically coupled. For G-type antiferromagnetism, 
the Fe spins are antiferromagnetically coupled both out-of- 
plane and in-plane. The shaded box indicate the unit cell of 
two formula unit supercell. 



Two sets of supercell calculations were carried out, one 
with two formula unit and another with eight formula 
unit. The two formula unit supercells with two inequiv- 
alent Fe atoms can accommodate the ferromagnetic spin 
alignment of Fe spins and the A type antiferromagnetic 
spin alignments of Fe spins. The eight formula unit su- 
percells with eight inequivalent Fe atoms in the unit cell, 
in addition to FM and A type AFM, can accommodate 
G type antiferromagnetic ordering of Fe spins (see Fig.l) 

For extraction of a few-band, tight-binding Hamilto- 
nian out of full DFT calculation which can be used as 
input to multi-orbital, low-energy Hamiltonian based cal- 
culations, we have carried out NMTO-downfolding calcu- 
lations. Starting from a full DFT calculations, NMTO- 
downfolding arrives at a few-orbital Hamiltonian by inte- 
grating out degrees which are not of interest. It does so 
by defining energy-selected, effective orbitals which serve 
as Wannier-like orbitals defining the few-orbital Hamilto- 
nian in the downfolded representation. NMTO technique 
which is not yet available in its self-consistent form re- 
lies on the self-consistent potential parameters obtained 
out of linear muffin-tine orbital (LMTO)^ calculations. 
The results were cross-checked among the plane wave and 
LMTO calculations in terms of total energy differences, 
density of states and band structures. 

The multi-orbital, low-energy Hamiltonian that is as- 
sumed to capture the essential physics of SFMO, should 
consist of the following ingredients: (i) a large core spin 
at the Fe site, (ii) strong coupling on the Fe site between 
the core spin and the itinerant electron, strongly prefer- 
ring one spin polarization of the itinerant electron, and 
(iii) derealization of the itinerant electron on the Fe-Mo 
network. 

From the above considerations, the representative 
Hamiltonian is given by: 
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-t F M fi*, a m 3<r,a ~ *MM ^ TO L,a m j<r,a 

<ij>(T,a <ij ><r,a 

X! fi<T.afjcr,a + Sj ■ f iot 3 aj3 fip (1) 

<2J>cr,Q 

The /'s refer to the Fe sites and the m's to the Mo 
sites. tpM, t m m j £ff represent the nearest neighbor Fe- 
Mo, second nearest neighbor Mo-Mo and Fe-Fe hoppings 
respectively, the largest hopping being given by t F M- & 
is the spin index and a is the orbital index that spans 
the t 2s manifold. The difference between the t2 9 levels of 
Fe and Mo, A = e Fe — ^Mo, defines the charge transfer 
energy. Since among the crystal-field split d levels of Fe 
and Mo, only the relevant t 2g orbitals are retained, the 
on-site and hopping matrices are of dimension 3x3. 
The Si are 'classical' (large S) core spins at the B site, 
coupled to the itinerant B electrons through a coupling 
J ^> t FM . 

Given the fact that J ^> tpM, the Hamiltonian of 
Eqn [l] can be cast into form appropriate for J — > oo. 
This gives the following Hamiltonian, with 'spinless' Fe 
conduction electrons and Mo electrons having both spin 
states. 

H = t FM ]T {sin{°^)f} a m^ a ~ e l ^cos{ e -^)f} a m na ) 

<ij>a 

+h.c. + t MM Y m \ocJ n i<ya 

<ij> 

+t FF ««(%/2)(/LiW) 

<ij> 

+ CFe fi<*fi<* + 6M ° m L a m iaa (2) 
i iaa 

This is the lowest energy Hamiltonian. There is no 
longer any 'large' coupling in the Hamiltonian, and the 
number of degrees of freedom has been reduced to three 
per Fe site and six per Mo, compared to original prob- 
lem with six degrees of freedom at both Fe and Mo sites. 
uiji and m,j-\ hop to different conduction electron projec- 
tions at the neighboring Fe sites so the effective hopping 
picks up a 9i,cj>i dependent modulation. For example, 
9 = 0, <j> = 0, corresponds to FM configuration with all 
Fe core spins being up. Since the spin S is large and can 
be considered classical, one can consider different spin 
configurations (ferro, antiferro and disordered) and di- 
agonalize the system in real space, to obtain variational 
estimates of the ground state, and its stability. 

III. RESULTS 

A. Total Energy, Electronic Structure and Relative 
Stability of Magnetic Phases 

Sr 2 FeMoOg crystallizes in body centered tetragonal 
space group with I4/mmm symmetry. The crystal struc- 



ture of SFMO is well characterized. The crystal structure 
of La-doped S^FeMo Og on the other hand is controver- 
sial. Some of the studjBEH reports that though I4/mmm 
symmetry is retained for small doping of La, for doping 
beyond x — 0.4 or so, the symmetry changes to P2i/n. 
The other measurements^ however reports that all com- 
pounds of Sr 2 _ x La a ,FeMo0 6 for x = 0, 0.25, 0.5 and 1.0 
crystallize in I4/mmm symmetry. Unfortunately, the in- 
formation of the detail crystal structure data are limited 
due to the facts that a) the compounds till now have 
been synthesized only for La concentrations less than or 
equal to 1 b) increasing concentration of La leads to in- 
creased disorder which prohibits accurate measurement 
of the underlying symmetry. While in the following, we 
have primarily reported the results assuming I4/mmm 
symmetry, we have also carried out calculation for P2i/n 
symmetry for the end member, La2FeMoOg (LFMO). 
The crystal structure corresponding to P2i/n symmetry 
for La2FeMo06 was obtained starting with initial param- 
eters of x — 0.4 as reported in ref.^and then performing 
total energy optimization of the initial structure. The 
P2!/n symmetry structure has been found to be energet- 
ically lower in energy by 90 meV than the correspond- 
ing I4/mmm symmetry structure. However, as described 
later, the primary conclusion of our results is found to 
remain unaffected by this possible change of symmetry. 
Table I. shows the theoretically optimized crystal struc- 
tures obtained using plane wave basid^O of SFMO, and 
that of LFMO assuming I4/mmm symmetry as well as 
P2!/n symmetry. 

The volume for LFMO is found to expand with re- 
spect tot hat of SFMO, in agreement with experimental 
increasing volume with increased La dop- 
ing. Assuming 14 / mmm symmetry, as is seen from Table 
I, the internal parameters corresponding to oxygen posi- 
tions, which are the only free parameters within I4/mmm 
space group, change little upon changing Sr by La. The 
unit cell volume for various intermediate members of the 
series obtained by interpolation from the optimized lat- 
tice parameters of the end members using Vegar's law, 
120.99 A 3 for SFMO, 123.27 A 3 for Sri. 5 La . 5 FeMoO 6 
and 125.56 A 3 for SrLaFeMo0 6 , agree well with the 
experimental data available for I4/mmm symmetry in 
terms of volume expansion ; given by 121.4 A 3 , 124.0 A 3 
and 124.88 A 3 respectively^. The crystal structure for 
the doped compounds in the assumed I4/mmm symme- 
try for the intermediate concentration values are, there- 
fore, obtained by using Vegard's law for interpolation of 
cell parameters keeping the atomic positions fixed. 

In the next step, we have carried out total energy cal- 
culations of Sr2_a;La x FeMo06 in I4/mmm symmetry for 
the FM alignment of Fe spins and the AFM alignment 
of Fe spins, which for a two formula unit supercell is 
of A type (see Fig. 1). The energy difference between 
FM and AFM-A spin configuration per formula unit as 
a function of La concentration is plotted in Fig. 2. Cal- 
culations have been carried out both within GGA and 
GGA+U. Focusing on to GGA results first, as is evi- 
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Table I: Optimized cell parameters and the atomic positions for S^FeMoOg and La2FeMo06- Fe and Mo ions are situated at 
the high symmetry Wykoff positions 2a and 2b, given by (0,0,0) and (0.0, 0.0, 0.5) respectively. For I4/mmm symmetry Sr/La 
also sites in the high symmetry Wykoff position given by (0.5, 0.0, 0.25) but sits in a general position for P2i/n symmetry. 



dent from Fig. 2, the stability of the FM phase with re- 
spect to AFM configuration is gradually reduced as the 
La concentration is increased. As the concentration is 
increased beyond x=1.5 or so, the FM phase becomes 
unstable and the AFM phase becomes the ground state, 
in agreement with prediction of model calculation J 27 H 26 | 
The total and magnetic moments at Fe and Mo sites, 
as obtained within GGA, are listed in Table II. The net 
magnetic moment at the FM phase reduces as the La 
concentration is increased, which is due to the increased 
moment at the Mo site. Such behavior has been also ob- 
served in experimenlP^. Especially, photoemission stud- 
ies have confirmed that electron injection occurs at the 
Mo site, increasing the moment on that site^l While 
the moment at the Fe site stays more or less the same 
between ferromagnetic and antiferromagnetic phase, the 
magnetic moment at the Mo site is found to be systemat- 
ically smaller in the AFM phase compared to FM phase. 

In order to check the influence of the possible change 
of crystal symmetry that may happen between SFMO 
and LFMO, we have also calculated the total energy dif- 
ference between FM and AFM-A spin configurations, as- 
suming LFMO in P2i/n symmetry with theoretically op- 
timized structure. The calculated Efm - ^afm-A came 
out to be 0.094 eV per formula unit, confirming the sta- 
bilization of AFM phase for LFMO. While the possible 
change of crystal symmetry from I4/mmm to P2i/n for 
La rich samples is expected to change the precise La con- 
centration at which FM to AFM transition happens, the 
general trend of AFM phase becoming progressively more 
favorable upon increasing La doping therefore would re- 
main hold good. 

Fig. 3 shows the GGA density of states correspond- 
ing to FM phase of SFMO, LFMO and the doped com- 
pounds, SrLaFeMoOg and Srg.sLa! 5 FeMoOg in I4/mmm 
symmetry. Focusing on the well-studied^ DOS of 
SFMO, we find that the Fe d states are nearly full 
(empty) in the majority (minority) spin channel while 
the Mo d states are nearly empty in the majority spin 
channel and partially filled in the minority spin channel. 
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Table II: Magnetic moments at Fe and Mo sites, and 
the total magnetic moment in FM and AFM-A phase 
of Sr2- a: La 2 ,FeMo06 in a two formula unit calculation. 
S3LFMO, SLFMO, SL3FMO refer to Sn.sLao.sFeMoOe, 
SrLaFeMoOe and Sro.sLai.sFeMoOe respectively. 



This is in conformity with the half metallic character of 
the compound and also with the nominal Fe 3+ and Mo 5+ 
valences. Due to the octahedral oxygen surrounding of 
Fe and Mo atoms, the Fe d and Mo d states are split 
up into t 2s and e g , the highly delocalized state cross- 
ing the Fermi level in the minority spin channel being of 
mixed Fe-t 2ff -Mo-t 2g character. The empty Mo-t 2g states 
in the majority spin channel is found to be highly local- 
ized giving rise to peaked structure positioned at about 
1 eV above the Fermi energy. As each of the Sr atoms is 
replaced by a La atom, one extra electron is introduced 
in the system which populates the hybridized Fe-t 2g -Mo- 
t 2g state in the minority spin channel, keeping the overall 
structure of the density of states intact. The Fermi level 
therefore progressively moves up like a rigid band fashion 
as x is increased and eventually hits the van Hove singu- 
larity of the Mo-t2 9 states in the majority spin channel. 
The FM solution becomes unstable at this point. This is 
schematically shown in the left panel of Fig. 5. Interest- 
ing the DOS corresponding to the mixed Fe-t 2g -Mo-t 2g 
character in the minority spin channel also exhibits the 
singularity at the same energy due to the essentially 2D- 
like nature of the hoppings between Mo-t 2g and Fe-t 2g 
Wannier functions as will be discussed in the following 
section. 
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Figure 2: The energy difference between FM and AFM-A 
phase plotted as a function of La concentration. The FM 
phase becomes unstable beyond a critical concentration of La 
both within GGA and GGA+U. 
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Figure 3: (Color online) The GGA density of 

states corresponding to SFMO, SrLaFeMoOe (SLFMO), 
Srn. 8 Lai. s FeMo0 6 (L3SFMO) and LFMO in ferromagnetic 
configuration. The density of states projected onto Fe, Mo 
and O are represented by solid black, green (grey) and shaded 
grey area. The upper and lower panels correspond to major- 
ity and minority spin channels. Zero of the energy is set at 
the GGA Fermi energy. 



Fig. 4 shows the density of states of SFMO, 
LFMO and the doped compounds, SrLaFeMoOg and 
Sr 5 La! 5 FeMo06 in the antiferromagnetic A phase, as 
calculated within GGA. In the two formula unit super- 
cells, there are two inequivalent Fe atoms, Fel and Fe2, 
whose spins are antiferromagnetically oriented. The ma- 
jority channel of Fel therefore is identical to the minor- 
ity channel of Fe2 and vice versa. The induced moments 
at two inequivalent Mo sites also become antiferromag- 
netically aligned, giving rise to a net AFM arrangement 
with a zero total moment. Shown in Fig 4, are therefore, 
the partial DOS corresponding to one of the sublattice 
since that of the other sublattice is identical with major- 
ity and minority spins reversed. We find that the Mo-Fe 
hybridized state crossing the Fermi level, has a three peak 
van Hove structure. This arises because of the fact that 
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Figure 4: (Color online) The density of states correspond- 
ing to SFMO, SrLaFeMoOe (SLFMO), Sro.sLai.sFeMoOe 
(L3SFMO) and LFMO in the A-type antiferromagnetic con- 
figuration. The density of states projected onto Fe, Mo and O 
are represented by solid black, green (grey) and shaded grey 
area. Zero of the energy is set at the GGA Fermi energy. 



due to creation of sublattices in the AFM phase, the Mo 
hopping becomes restricted to a reduced dimension as 
the Mo electrons can effectively hop to Fe sites with a 
specific orientation of Fe spins and not in another. In- 
terestingly, such a three peak structure formation is also 
seen in model calculation (see Fig-3 of ref[28]). As found 
in the case of FM DOS, the gross features of the density 
of states remain unchanged with the La doping apart 
from the upward shift of the Fermi energy. Reaching 
LFMO, the Fermi level lands up in the dip of the three 
peak structured DOS, justifying the stability of the anti- 
ferromagnetic phase, as shown in the schematic diagram 
of Fig. 5. 

The antiferromagnetic state becomes energetically fa- 
vorable, when the filling is such that it starts populating 
the Mo states in the majority spin channel of the FM 
DOS, which is highly localized due to the strong prefer- 
ence of the Mo-Fe hopping in one spin channel and not 
in another. The antiferromagnetic configuration of Fe 
spins, on the other hand, allows both Mo down spin as 
well as up spin electron to hop, albeit in different sublat- 
tices, thereby stabilizing the AFM phase through kinetic 
energy gain. 

In order to check the influence of the missing correla- 
tion effect in GGA, we have also carried out GGA+U cal- 
culations with a typical U valued of 4 eV and J value of 
1 eV, applied at the Fe site. The calculated energy differ- 
ence between FM and AFM-A configuration as a function 
of La doping is shown in Fig 2, along with GGA results. 
The application of U is found to increase the relative sta- 
bility of AFM phase due to the increased superexchange 
contribution to antiferromagnetism in addition to kinetic 
energy driven antiferromagnetism. 

In Fig [6] we show the GGA+U DOS for LFMO, plot- 
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Figure 5: (Color online) Schematic diagram showing the 
mechanism stabilizing the AFM phase over FM phase. As 
the La doping is increased, the Fermi level (E/) shifts towards 
right. 

ted for both FM and AFM-A phases. It is observed that 
the gross features of the DOS close to Fermi energy, re- 
main similar to GGA: in particular, the Fermi energy in 
the FM phase remains pinned to the unoccupied Mo t^g 
level in the FM phase. However, the hybridization be- 
tween the Fe and Mo decreases. Nevertheless, the anti- 
ferromagnetic state is still found to have a finite density 
of states at Fermi energy, signifying the dominance of 
kinetic energy driven contribution over that of superex- 
change. 

As already mentioned, considering the two formula 
unit supercell, the possible AFM arrangement that can 
be achieved is of A type. In order to achieve the G type 
AFM ordering involving both in-plane and out-of-plane 
AFM ordering, one needs to increase the size of the su- 
percell to at least eight formula unit. Eight formula unit 
supercells also allow to probe the concentration range 
intermediate to x=1.5 and x=2.0, the region where the 
crossover from FM to AFM happens. Since the quali- 
tative conclusions remain unchanged between GGA and 
GGA+U, the eight formula unit calculations were car- 
ried out only for GGA. The energy differences per for- 
mula unit obtained for different concentrations of La be- 
tween FM and AFM-A, and between FM and AFM-G 
phases are listed in Table III. As found in the calculations 
with two formula unit, the stability of the FM phase is 
found to gradually decrease as the La concentration in- 
creases. Among the two antiferromagnetic phases, the 
G-type AFM is found to be energetically very close to A 
type AFM phase, with G type AFM being the ground 
state at the end limit of doping i.e. for LFMO. 

B. Determination of low-energy, few orbital 
Hamiltonian by NMTO-downfolding 

In order to probe the variation of La concentration in 
a continuous manner, it is perhaps more convenient to 
adopt a low-energy Hamiltonian approach. This would 
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Figure 6: (Color online)DOS for LFMO in FM and AFM-A 
phase, using GGA+U 
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Table III: Total energy differences per unit formula in eV be- 
tween FM and AFM-A, and between FM and AFM-G for 
various doping of La, as obtained within eight formula unit 
supercell calculations. 



also allow one to calculate the physical properties like 
magnetic transition temperatures, transport, spin wave 
spectra in a much more manageable way. For this pur- 
pose, however, it is essential to construct a realistic, low- 
energy Hamiltonian. We have used for this purpose, the 
NMTO-downfolding technique. For the present problem, 
we have derived a Fe-t2 S - Mo-t2 9 only Hamiltonian by 
integrating out all the degrees of freedom other than Fe- 
t<i g and Mo-t 2j . Calculations were carried out both in 
the spin-polarized and non spin-polarized form. First of 
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all, Fig.[7]illustrates the driving mechanism of magnetism 
in this class of compounds^. The top panels show the 
on-site energies of the real-space Hamiltonian defined in 
downfolded effective Fe-Mo basis for SFMO and LFMO 
in a spin polarized calculation. As is seen, the t2 S levels 
of Mo appear in between the exchange split Fe d states. 
Upon switching on the hybridization between Fe-d and 
Mo-t2 ff , states of same symmetry and spin interact. As 
a result, Mo-t2 9 up spin states are pushed up in energy 
and Mo-t2 S down spin states are pushed down in energy, 
introducing a renormalized, negative spin splitting at the 
Mo site. The normalized spin splitting at Mo site is esti- 
mated by massive downfolding procedure by keeping only 
Mo-t2 ff states active in the basis, as shown in the right 
half on the top panels in Fig. [7] . We note that this to be 
true for both SFMO and LFMO. This in turn, once again, 
reconfirms the hybridization driven mechanism to be op- 
erative both in SFMO and LFMO, the only difference 
being in the carrier concentration. This is in contrast to 
Sr2FeMo06 where W t2g-0 hybridized levels are pushed 
above the exchange split Fe d levels. The increase in the 
number of conduction electrons for LFMO compared to 
SFMO, is reflected in the spin splitting at Mo site before 
switching of the hybridization, to be about three times 
larger in LFMO (0.37 eV) compared to that of SFMO 
(0.13 eV). The bottom panels of Fig. |7]show the plots of 
Wannier functions of the massively downfolded Mo-t2 S in 
the down spin channel which demonstrates the hybridiza- 
tion between Mo-t2 9 and Fe-t2 S states. 

Table IV. shows the hopping interactions between Fe 
and Mo, obtained in the basis of Fe and Mo t2 9 Wannier 
functions constructed by NMTO-downfolding technique. 
The numbers inside the bracket are that of LFMO while 
those outside are that of SFMO. The examination of the 
hopping table reveals that the nearest neighbor Fe-Mo 
hopping to be strongest, as expected. The second nearest 
neighbor Mo-Mo hopping is half as strong as the near- 
est neighbor Fe-Mo hopping, while the second nearest 
neighbor Fe-Fe hopping is about one fifth of the near- 
est neighbor Fe-Mo hopping. The out-of-plane hoppings 
which are of dd5 kind are order of magnitude smaller 
than the in-plane dd7r kind of hopping, while inter-orbital 
hoppings are found to be negligibly small (less than 0.01 
eV). This makes the hopping essentially two dimensional, 
as commented earlier. As is also evident, by replacing 
Sr by La, the essential material specific parameters of 
the low-energy Hamiltonian, as given in Eqn. 1 changes 
very little. This is shown pictorially in Fig [7] where it is 
found that the relative energy positions of the i 2 gj. levels 
of Fe and Mo change very little in going from SFMO to 
LFMO. In the solution of low-energy Hamiltonians, to 
be described in the next section, calculations are there- 
fore carried out assuming the hopping parameters cor- 
responding to SFMO and varying the carrier concentra- 
tion in a rigid-band fashion. The charge transfer energy 
between Fe-t2 ff and Mo-t2 9 has been found to differ by 
about 5% which has been taken into account in these 
calculations. 



SFMO LFMO 




Figure 7: (Color online) Top panels: Positioning of vari- 
ous energy levels as obtained by NMTO- downfolding calcu- 
lation before and after switching on the hybridization be- 
tween the magnetic and nonmagnetic ions. Bottom panels: 
Effective Mo-t2 9 Wannier orbitals corresponding to massively 
downfolded NMTO Hamiltonian in the down spin channel. 
Shown are the orbital shapes (constant-amplitude surfaces) 
with lobes of opposite signs colored as blue (dark grey) and 
cyan (light grey). The central part of the Wannier orbitals 
are shaped according to Mo-t2 S symmetry, while the tails are 
shaped according to Fe-t2 9 and O-p symmetries. Significant 
amount of weights are seen in O and Fe site which reflects 
the strong hybridization between Fe, Mn and O. For LFMO, 
finite weights are seen also at La sites, occupying the hollows 
formed between Mo-0 and Fe-0 bonds, which is of La 3z 2 
character. 



C. Calculations of Magnetic Phase Diagram and 
Magnetic Transition Temperatures in terms of 
Low-energy Hamiltonian 

The exact diagonalization of the low-energy Hamilto- 
nian, as given in Eqn. 2 has been carried out for finite 
size lattice of dimensions 4x4x4, 6x6x6 and 8 x 
8x8. The hopping parameters and the onsite energies 
were taken out of DFT calculations, as listed in Table 
IV. For convenience of calculation, we have neglected the 
small tetragonality reflected in the parameters listed in 
Table IV. The dominant hopping interaction which is be- 
tween nearest-neighbor Fe and Mo is found to be of the 
order of 0.3 eV, while the spin exchange splitting at Fe 
site as shown in Fig. [7] is of order of the 3 eV, an order of 
magnitude larger than the dominant hopping interaction. 
This justifies the assumption of J — > oo limit as adopted 
in Eqn. 2. This makes the rank of the Hamiltonian to 
be diagonalized as 9/2 x N 3 for a N x N x N lattice^. 

The energy difference between ferromagnetic configu- 
ration and G-type and A-type antiferromagnetic config- 
uration of Fe spins as a function of carrier concentra- 
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Direction 


xy,xy 




xz,xz 


INN 
(Fe-Mo) 


[100] 
[010] 
[001] 


-0.26 (-0.26) 
-0.26 (-0.26) 
-0.02 (-0.04) 


-0.02 (-0.04) 
-0.26 (-0.26) 
-0.26 (-0.25) 


-0.26 (-0.26) 
-0.02 (-0.04) 
-0.26 (-0.25) 


2 NN 
(Fe-Fe) 


[110] 
[101] 
[011] 


-0.05 (-0.06) 
0.00 (0.00) 
0.00 (0.00) 


0.01 (0.00) 
0.01 (0.00) 
-0.04 (-0.06) 


0.01 (0.00) 
-0.04 (-0.06) 
0.01 (0.00) 


2 NN 
(Mo-Mo) 


[110] 
[101] 
[011] 


-0.11 (-0.12) 
0.01 (0.01) 
0.01 (0.01) 


0.00 (0.01) 
0.00 (0.01) 
-0.11 (-0.12) 


0.00 (0.01) 
-0.11 (-0.12) 
0.00 (0.01) 


3NN 
(Fe-Mo) 


[111] 


0.01 (0.00) 


0.00 (0.00) 


0.00 (0.00) 


4 NN 
(Fe-Fe) 


[100] 
[010] 
[001] 


0.01 (0.01) 
0.01 (0.01) 
0.01 (0.00) 


0.01 (0.00) 
0.01 (0.01) 
0.01 (0.01) 


0.01 (0.01) 
0.01 (0.00) 
0.01 (0.01) 


4 NN 
(Mo-Mo) 


[100] 
[010] 
[001] 


0.01 (0.03) 
0.01 (0.03) 
0.01 (0.00) 


0.01 (0.00) 
0.01 (0.03) 
0.01 (0.03) 


0.01 (0.03) 
0.01 (0.00) 
0.01 (0.03) 


5 NN 
(Fe-Mo) 


[110] 
[101] 
[011] 


-0.01 (-0.01) 
0.00 (0.01) 
0.00 (0.00) 


0.00 (0.01) 
0.00 (0.00) 
-0.01 (-0.01) 


0.00 (0.00) 
-0.01 (-0.01) 
0.00 (0.01) 



Table IV: Hopping matrix elements in eV between Fe-t2 ff and 
Mo-t2 9 . Only the hopping matrix elements of magnitude 
larger than 0.01 eV are listed. The onsite matrix elements are 
given by 0.005 (0.008) eV, 0.0 (0.0) eV, 0.0 (0.0) eV for Fe- 
xy, Fe-j/z and ¥e-xz respectively, and 1.018 (1.057) eV, 1.007 
(1.053) eV, 1.007 (1.053) eV for Mo-xy, Mo-yz and Mo-iz re- 
spectively. All numbers inside the bracket are for LFMO and 
those outside are for SFMO. The energies for a given com- 
pound is measured with respect to the lowest energy state. 
The small differences between numbers involving xy and that 
of yz and xz reflect the tetragonality present in the systems. 

tion is plotted in Fig. [8] . The negative values of the 
energy differences indicate the stability of the ferromag- 
netic phase while the positive values indicate the stabil- 
ity of the antiferromagnetic phase. The cross-over hap- 
pens for a value of conduction electrons equal to about 
~ 2.6, corresponding to x = 1.6, which agrees well with 
the results of 8 formula unit supercell calculations, given 
the assumption of infinite Hund's coupling at Fe site and 
the finite size effect. This agreement is nontrivial, since 
the effective Hamiltoanian has only 12 spin-orbitals, and 
hence 12 bands, as compared to the 500 band calcula- 
tion with 8-formula unit supercells. This in turn, vali- 
dates the construction of low-energy model Hamiltonian 
as given in Eqn. (2), in terms of correct identification of 
the essential contributing terms. This gives us confidence 
in the constructed low energy model Hamiltonian, which 
can henceforth be used to calculate many other proper- 
ties like conductivity, susceptibility, magnetoresistance, 
including at finite temperature, which are not easily ac- 
cessible within DFT. 

As an example, we have used the solutions of the low- 
energy Hamiltonian to calculate the magnetic transition 




No of conduction electron 



Figure 8: The energy differences between the FM and G- 
type AFM phase (solid line) and the FM and A-type AFM 
phase (dashed line) plotted as a function of the number of 
conduction electrons, as obtained by exact diagonalization of 
the low-energy Hamiltonian for a 8 x 8 x 8 lattice. Only the 
region outside the hashed regions, from carrier concentration 
1 to 3 is of relevance for S^-aTa^FeMoOe- 




Number of conduction electrons 



Figure 9: The ferromagnetic T c (solid line) and the antiferro- 
magnetic transition temperature Tjv (dashed line) plotted as 
a function of the number of conduction electrons, as obtained 
by exact diagonalization of the low-energy Hamiltonian for a 
8x8x8 lattice. As in Fig pi, the region outside the hashed 
regions, from carrier concentration 1 to 3 is of relevance for 
Sr 2 _^La :c FeMo06. 



temperatures by calculating the difference between the 
paramagnetic phase and the relevant magnetic phases. 
The paramagnetic phase was simulated as disordered lo- 
cal moment calculations, where the calculations were car- 
ried out for several ( ~ 50) disordered configurations of Fe 
spin and were averaged to get the energy corresponding 
to paramagnetic phase. We note that, such a calcula- 
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tion would have been rendered extremely difficult within 
ab-initio owing to the computational time involved using 
large supercells, and also averaging them over myriad 
configurations. Fig. [9] shows the transition temperatures 
plotted as a function of carrier concentration. The ferro- 
magnetic transition temperature at carrier concentration 
of 1, which corresponds to SFMO compound, is found to 
be 360 K in comparison to measured value of 410 KP^. 
The ferromagnetic T c is found to decrease upon increas- 
ing La concentration, and finally becomes zero. Upon 
suppression of ferromagnetic T c , the transition tempera- 
ture of the antiferromagnetic phase, T^ starts growing, 
hitting a maximum value for the end member, LFMO. 



IV. SUMMARY AND OUTLOOK 

Using the combination of first-principles DFT calcula- 
tions and exact diagonalization calculations of low-energy 
Hamiltonians, we showed that the La doped Sr 2 FeMo06 
compounds become progressively more unstable towards 
ferromagnetism as the La concentration is increased. For 
the La rich members of Sr 2 _ a; La :I .FeMo06 series with 
x > 1.6, the ground state becomes antiferromagnetic. 
This antiferromagnetic phase is found to be governed 
by the kinetic energy driven mechanism as operative 
in SFMO and achieved by change in carrier concentra- 
tion. In contrast to the super-exchange driven antifer- 
romagnetic phase found in case of double perovskite like 
SraFeWOePH this antiferromagnetic phase presumably is 
metallic. Our DFT calculations found antiferromagnetic 
solutions with finite density of states at Fermi energy. 
The preliminary calculations of the matrix elements of 
the current operator also turned out to be non-zero. This 
will be taken up with more rigor in near future. 

It is interesting to compare our results with Dynamical 
Mean Field Theory (DMFT) calculations done by Chat- 
topadhyay and Millis^, using a one-band model Hamil- 
tonian. This was, however, a single site calculation, and 
hence there was no possibility of capturing an antiferro- 
magnetic phase. Their T c vs N plot for the ferromag- 
netic phase, however, was very similar to ours, as shown 
in Fig. 9, in the sense that the T c first increased, and 
then decreased with increasing filling, finally, becoming 
at a filling close to 3. The additional and the most im- 
portant finding of our study is that our calculations also 
demonstrate the cause of the vanishing T c : namely the 
emergence of the AFM phase. 

There are however, several important issues which 
needs to be considered. Formation of pure La2FeMoOg 
to best of our knowledge is not been reported in liter- 
ature, which presumably is due to relative scarcity of 
Mo 3+ ions in octahedral environment. La rich SFMO 
samples, therefore seem more promising candidates for 
exploration of the antiferromagnetic phases. 

Another important issue to bother about is the issue 
of antisite disorder, which has not been considered in 
our study. The necessary conditions of formation of dou- 



ble perovskites with ordered, rock salt arrangement of 
B and B' transition metal ions are that the size differ- 
ence between B and B' ions should be sufficiently large 
as well the nominal charge difference. With the increase 
of La concentration, the extra doping electrons populate 
the Mo t2 g down spin sub-band crossing the Fermi en- 
ergy. As a result, the Mo +5 nominal valence in SFMO 
becomes Mo +5_m in the doped compounds, m being the 
number of doped electron with a maximum value of 2. 
This decreases the charge difference between Fe +3 and 
Mo. This is expected to be detrimental to the ordering, 
though the ionic radii difference between Fe 3+ ( 0.645 A) 
and Mo 3+ (0.69 A) is larger than that between Fe 3+ and 
Mo 5 + (0.65 A). The study on SFMO in this contexfP, 
find that even for a disordered sample, as given by X-ray 
study, the local ordering is maintained with a domain 
structure. Annealing conditions can give rise to domain 
structures with varying sizes of the domain. This gives 
us the hope in the observation of the antiferromagnetic 
phase in the La rich SFMO samples. Attempts are al- 
ready being made to prepare these overdoped samples 
locally, and preliminary results suggest existence of mag- 
netic phases different from ferromagnetic phased On 
theoretical front, effect of disorder has been studied by 
Alonso et. alP^. Within a variational mean field frame- 
work, they found that the filling/doping at which the Tc 
goes to zero increases upon increasing antisite disorder. 
This means that antisite disorder stabilizes the ferromag- 
netic phase. In other words, the antisite disorder is ex- 
pected to reduce the stability of the antiferromagnetic 
phase, which is also seen in our preliminary calculations. 
We wish to carry out systematic study of the antisite 
disorder in future, keeping in mind possible formation of 
domain structures. 

Finally, within the kinetic energy driven mechanism, 
the ferromagnetism gets destabilized and the antiferro- 
magnetism wins when the carrier concentration reaches 
such a value that the B' d states gets filled up in one 
spin channel and tries to populate the other spin chan- 
nel. Such situation is encountered also in case of an- 
other double perovskite, Sr 2 CrOs06- Os being in nomi- 
nal 5+ state, is in d 3 configuration with completely full 
Os t 2ff states in the down spin channel and lies within ex- 
change split energy levels of Cr-t2cp^, a case very similar 
to La2FeMoOg. Sr 2 CrOs06, in contrast to above expec- 
tation, however stabilizes in ferromagnetic configuration 
of Cr spins. In this context, we found that the delicate 
balance between FM and AFM, is governed by the extent 
of hybridization between the localized B site and delocal- 
ized B' site. For S^CrOsOe, due to the movement of the 
Os t2 9 within the exchange split energy window of Cr- 
t2 9 in comparison to that of Mo t2 9 within the exchange 
split energy window of Fe-d, the hybridization effect is 
weakened and also the finite spin-orbit (SO) coupling at 
Os site mixes the up and down spin channels, causing 
possibly the energy gain due to antiferromagnetism to 
be reduced^. 

We believe our study will stimulate further experimen- 
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